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June  1980  Updates  to  SOQ80128 


INTRODUCTION 

This  document  defines  the  changes  made  to  the  SOQ  code  (SOQ80128) 
between  January  and  June  of  1980.  The  changes  either  correct  short¬ 
comings  found  in ~ the  code  or,  more  usually,  document  the  increased 
capability  being  continually  built  into  the  code.  The  SOQ  code  is 
maintained  as  SOQ80128  June  PL, ID  *  AFLOJRA  as  a  NOS/BE-1  CDC  update 
format  file. 

UPDATES 

1 .  *ID  FIXZRN 

This  update  redefines  the  coeffiencets  to  be  input  to  the  Zernike 
subroutine.  This  new  convention  is  more  physically  meaningful  in  that, 
at  least  for  lower  orders,  the  coefficients  are  in  waves.  For  example, 
to  impose  one  wave  peak  to  peak  of  defocus  (P^)  on  a  beam,  one  would 
input  P(4)-l.  The  phase  applied  is  now: 

♦  d.J)  =  |Pk*zk(i,j) 

The  subroutine  affected  is  ZERN.  This  update  does  not  effect  the  rest 
of  the  code. 

2.  *ID  FIXJTR 

This  update  ensures  a  correct  definition  of  DF  in  subroutine  JITRBG 
since  when  JITRBG  is  called  from  subroutine  QUAL,  the  X-coordinate 
array  contains  RX/D  coordinates,  not  the  spatial  coordinates. 

Only  one  line  of  the  code  is  affected  by  this  update. 

3.  *ID  ROTZRN 

Due  to  different  coordinate  system  orientations  for  data,  it  became 
necessary  to  allow  for  this  variation  within  subroutine  ZERN. 

Define  the  data  x  and  y  coordinates  to  be  XROT  and  YROT,  and  the  SOQ 
x  and  y  coordinates  to  be  XIN  and  YIN.  The  rotation  angle  is  then 
defined  to  be  6  (in  radians). 
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COS ROT  *  C0S(6) 

SINROT  =  SIN(0) 

XROT  =  XIN  x  COS ROT  +  YIN  x  SINROT 
YROT  =  -XIN  x  SINROT  +  YIN  x  COSROT 

Application  of  Zernike  polynomials  to  and  SOQ  point  located  at 
(XIN,  YIN)  would  then  be  calculated  using  Z(XR0T,  YROT),  The  possi¬ 
bility  of  axis  flips  are  also  accounted  for  and  are  flagged  by 
FLIPX  or  FLIPY  not  equal  to  zero.  Namelist  ZERNS  is  modified  to 
include  FLIPX,  FLIPY  and  the  rotation  angle  (in  degrees)  ZTHETA.  No 
common  was  modified.  This  update  modified  only  subroutines  GDL  and 
ZERN. 


I 


rcFc** •  ic=  aPLCSSi  ♦  scceci2e»ic=ARLCvRA 

ICE  NT  FI  >2RN 


f  t  ge  i 


♦  /  2  ern 
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C  PLIFY  =  1.  RESLLTS  IN  A  FLIF  APCLT  TFE  Y  AXIS  CPFCFE 
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SECTION  I 
INTRODUCTION 

As  a  prerequisite  to  the  design  and  optical  integration  of  high  power 
laser  systems,  sophisticated  analytical  tools  are  needed  to  predict  overall 
system  performance  and  to  perform  design  tradeoff  studies.  Without  analytical 
guidance,  the  search  for  the  highest  level  of  safe  system  operation  becomes 
an  empirical  trial -and-error  process  involving  not  only  expensive  and  time- 
consuming  hardware  changes,  but  also  risk  to  critical  system  components. 

One  useful  analytical  tool  which  has  been  developed  for  the  design  and 
integration  of  electro-optical  systems  is  the  system  optical  quality  (SOQ) 
computer  code.  This  code  is  designed  to  guide  the  optical  designer  in  the 
comparison  and  assessment  of  system  optical  performance  by  providing  a 
physical  optics  model  of  the  system  which  traces  the  beam  from  its  point  of 
origin  in  the  resonator  through  the  optical  train  into  the  far  field. 

To  predict  overall  system  performance,  detailed  knowledge  of  the  optical 
field  (phase  and  amplitude  distributions)  is  needed  at  every  axial  station  in 
the  optical  train.  A  wave-optics  theory  is  mandatory  since  the  resonator  mode 
structure,  aperture  losses,  optical  feedback,  flux  spillage,  etc.,  are 
governed  principally  by  diffraction.  Moreover,  system  components  cannot  be 
analyzed  independently  since  the  beam  distortions  induced  by  each  component 
depend  on  the  shape  and  magnitude  of  the  incident  field  distribution. 

The  purpose  of  this  report  is  to  provide  a  fundamental  description  of 
the  SOQ  code  structure  and  a  description  of  the  input  of  the  code  as  aoplied 
to  gas  dynamic  laser  systems.  This  description  allows  the  optical  designer  to 
access  and  use  the  code  to  model  a  wide  variety  of  resonator  and  pointer/ 
tracker  concepts.  This  report  does  not  attempt  to  describe  how  the  system 
analyst  should  lay  out  his  optical  design  problem.  Familiarity  with  the 
optical  components  and  with  the  concepts  of  scalar  wave  theory  and  funda¬ 
mental  optics  are  assumed.  Also  assumed  is  a  working  knowledge  of  Fortran  IV 
and  the  CDC  NOS/BE  operating  system  commands.  Accurate  use  of  the  SOQ  code 
requires  familiarity  with  References  1  through  8. 

This  report  is  divided  into  four  main  sections.  Section  II  describes  the 
general  structure  of  the  SOQ  code  and  establishes  a  correlation  between  the 
usual  optical  elements  encountered  in  the  optical  train/gas  dynamic  laser 
resonator  and  the  appropriate  SOQ  models.  Section  III  acauaints  the  user  with 
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the  individual  SOQ  subroutines  and  their  analytical  formulations  as  mani¬ 
fested  in  Fortran  within  the  SOQ  framework.  This  section  also  delineates  the 
inputs  required  to  excerise  the  subroutines.  Section  IV  familiarizes  the  user 
with  the  operation  of  the  SOQ  model.  This  section  contains  working  input 
modules  which  carry  the  user  through  the  usual  calculations  of  the  SOO  code 
from  input  generation  to  loaded  cavity  calculations. 

An  additional  section,  Section  V,  has  been  included  to  show  the  user  how 
to  access  and  modify  the  code. 
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SECTION  II 
CODE  DESCRIPTION 


1.  GENERAL  CODE  OVERVIEW 

Section  II  provides  a  guide  to  the  region  of  applicability  of  the  SOQ 
code  and  a  synopsis  of  the  overall  structure  and  accessibility  procedures  for 
the  code.  It  is  to  acquaint  the  user  with  the  way  the  computer  code  operates 
and  with  the  subroutines  supplied  with  the  code. 

It  is  assumed  that  the  user  has  unfolded  his  optical  system  into  an 
equivalent  optical  train  in  which  each  wavefront /medium  interaction  may  be 
easily  identified  (Ref.  1).  The  following  illustrates  the  concept:  Figure  1 
is  a  simplified  diagram  of  a  telescope.  Figure  2  is  a  representation  of  the 
same  telescope  expressed  as  an  equivalent  optical  train.  The  SOQ  user  must 
correlate  the  optical  elements  in  the  code  (the  equivalent  system)  with  his 
optical  system. 


Figure  1.  Simplified  diagram  of  a  telescope. 


1.  Sziklas,  E.A.,  et  al.,  System  Optical  Quality  Study  Phase  I  --  Problem 
Definition.  Pratt  &  Whitney  Aircraft,  AFWL-TR,-?3-231 ,  Air  Force  Weapons 
Laboratory,  Kirtland  Air  Force  Base,  MM,  June  1974. 


Reference  Mirror- 


Figure  2.  Equivalent  optical  train  of  a  telescope. 

In  the  eauivalent  system,  the  actual  ODtical  train  is  replaced  by  an 
equivalent  optical  transmission  line  which  does  not  transmit  the  actual  beam, 
but  an  effective  unidirectional  beam.  In  the  wave  propagation  function  the 
propagation  algorithm  replaces  the  actual  beam  with  the  effective  unidirec¬ 
tional  beam  and  propagates  this  beam  through  the  entire  equivalent  optical 
train. 

The  optical  train  of  the  laser  system  under  investigation  is  divided 
into  a  series  of  linked  subsystems  A,  B,  C  .  .  .  .  etc.,  arranged  in  the 
sequence  encountered  by  the  beam  (Figure  3),  Laser  systems  employing  feed¬ 
back  (e.g.,  oscillators)  are  viewed  as  a  repetitive  chain  of  subsystems  in 
which  the  period  constitutes  one  round  trip  pass  through  the  resonator . 


Figure  5.  Optical  train  of  the  laser  system. 


Each  subsystem  A,  B,  C  .  .  .is  described  by  one  or  more  library  subrou¬ 
tines.  Subroutines  define  a  particular  optical  geometry  in  which  all  effects 
on  the  beam,  including  that  of  a  transmitting  inhomogeneous  medium,  are 
lumped  at  one  or  more  discrete  planes  along  the  optical  axis.  Thus,  each 
subroutine  defines  the  position  of  a  number  of  discrete  optical  elements 
(e.g.,  mirrors,  apertures,  gain/phase  segments)  along  the  optical  axis.  Free- 
space  propagation  is  prescribed  between  the  optical  elements.  Mirrors  and 
apertures  are  described  by  a  single  optical  element,  while  a  transmitting 
medium  is  described  generally  by  several  optical  elements.  The  number  of 
optical  elements  needed  is  governed  by  the  requirements  that  the  continuous 
amplitude  and  phase  changes  across  each  element  in  the  distributed  "loss” 
system  be  replaced  by  a  series  of  lumped  "losses'*  whose  amplitude  and  phase 
changes  are  small  relative  to  unity  and  2t,  respectively. 

The  input/output  format  of  each  subroutine  is  identical,  viz.,  an  N  x  N 
array  of  complex  numbers  representing  the  transverse  phase-amplitude  distri¬ 
bution  of  the  input  and  output  optical  fields.  In  this  discussion,  \T  is  an 
integer  defining  the  number  of  grid  points  used  to  specify  the  field  in  one 
transverse  dimension. 

Each  optical  element  is  characterized  by  a  complex  transmission  func¬ 
tion.  Thus,  if  u  (I,J)  denotes  the  field  incident  on  a  particular  element 
and  u+  (I,J)  the  output  field,  then 

u+(I,J)  »  u_(I,J)  (1) 

where  the  transmission  function  t  which  characterizes  this  optical  element 
can  be  written  in  the  general  form 

t(I,J)  *  A C I > J)  exp  |iA<fr(I,J)J  (2) 

Here,  A(I,J)  represents  the  amplitude  change  and  Ab(I,J)  the  phase  distortion 
induced  by  the  optical  element.  The  integers  I,J  =  1,  2,  .  .  .  N,  define  the 
two  transverse  grid  coordinates.  In  general,  a  subroutine  is  constructed  to 
define  each  transmission  function  and  then  stored  in  the  subroutine  library. 

In  the  wave  propagation  function,  an  effective  unidirectional  collimated 
beam  is  generated  and  propagated  along  the  equivalent  optical  transmission 
line.  The  beam  is  effectively  unidirectional  and  collimated  since  linear  and 
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quadratic  phase  components  (tilt  and  defocus)  are  normally  analytically 
removed  from  the  description.  Care  must  be  taken  to  properly  account  for 
the  apparent  transverse  displacement  of  optical  elements  when  tilt  is  analyt¬ 
ically  removed  from  the  wave.  The  resultant  collection  of  displaced  optical 
elements  defines  the  equivalent  optical  transmission  line,  as  illustrated  in 
Figure  4. 


Figure  4.  Equivalent  optical  transmission  line. 

The  system  performance  (e.g.,  beam  quality  at  various  specified  axial 
stations)  is  assessed  in  a  performance  evaluation  function.  This  function 
also  includes  initiation  of  any  desired  parametric  study  (e.g.,  cost/perform¬ 
ance  tradeoff  assessment)  in  which  the  wave  propagation  analysis  is  repeated 
for  each  specified  value  of  the  parameter. 

2.  REGION  OF  APPLICABILITY  OF  CODE 

The  SOQ  code  documented  in  this  report  is  applicable  to  calculations 
which  are  resolved  by  the  128  x  128  point  transverse  description  of  the  com¬ 
plex  field  behavior.  The  code  is  currently  outfitted  with  those  subroutines 
which  describe  a  flowing,  high  Mach  number  CO^-N-,  gas  dynamic  laser  system 
and  the  appropriate  optical  elements,  such  as  mirrors  and  apertures.  A  synop¬ 
sis  of  the  available  models  is  included  in  Section  II  of  this  report  and  a 
full  disclosure  is  included  in  Section  III.  The  128  x  128  point  code  is  well 
suited  to  the  description  of  problems  characterized  by  collimated  Fresnel 
numbers  15.  Use  of  the  128  x  128  point  code  for  propagation  at  higher 
Fresnel  numbers  has  been  shown  to  convolve  away  the  higher  spatial  frequency 
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information,  while  retaining  the  general  character  of  the  solution.  Thus,  the 
code  can  be  used  outside  the  Fresnel  number  region  described  above  without 
totally  erroneous  results.  Specific  requirements  are  set  in  the  following 
section  which  state  the  analytical  relationship  between  beam  geometric  site, 
required  calculation  region  site,  and  the  number  of  sampling  points. 

Though  currently  limited  to  123  x  123  points,  the  code  can  be  expanded 
to  256  x  256  points  for  loaded  cavity  calculations  on  the  Air  Force  Weapons 
Laboratory  (AFWL's)  current  Control  Data  Corporation  (CDC)  Cyber  176  com¬ 
puter,  The  SOQ  code  currently  uses  about  225gK  of  small  core  memory  (SCM)  and 
420gK  of  large  capacity  memory  (LCM)  to  load  and  execute.  The  basic  number  of 
point  requirements  is  specified  in  powers  of  2  to  achieve  the  most  efficient 
use  of  Fast  Fourier  Transform  (FFT)  techniques,  although  the  FFT  algorithm 
employed  for  the  propagation  algorithm  may  be  executed  with  any  integer 
radix,  but  with  some  loss  in  efficiency. 

The  SOQ  computer  code  installed  at  AFWL  includes  many  optical  transmis¬ 
sion  functions.  However,  the  basic  structure  of  the  code  allows  the  user  to 
include  whatever  additional  optical  elements  desired  by  appropriately  updat¬ 
ing  the  code.  A  detailed  description  of  subroutine  addition  is  included  in 
Section  V. 

The  computer  program  can  now  simulate  the  following  numbers  of  optical 
elements.  The  code  structure  limits  the  user  to  20  different  optical  ele¬ 
ments,  each  of  which  can  have  up  to  14  physical  parameter  input  arguments. 

Each  of  these  input  arguments  can  be  called  up  to  9  times  within  a  system 
setup.  Up  to  98  optical  elements  may  be  included  within  one  iterative  resona¬ 
tor  calculation  with  a  maximum  of  2  laser  flow  cavities  containing  5  gain- 
nhase  segments  per  cavity.  Each  of  these  gain-phase  segments  represents  a 
complex  transmission  function  which  couples  the  cavity  kinetics  and  flow 
field  inhomogeneities  with  the  propagating  optical  field.  Calculations  out¬ 
side  an  iterative  loop  limit  the  user  to  9  calls  to  the  same  subroutine.  All 
of  the  above  constraints  may  be  changed  by  the  user,  if  he  so  desires,  by 
redimensioning  the  deck.  Presently  the  program  also  provides  an  automatic 
procedure  for  storing  information  from  previous  calculations  and  restarting 
subsequent  calculations  based  on  these  results. 


3.  SAILING  REQUIREMENTS 

The  fundamental  SOQ  propagation  procedure  uses  a  numerically  efficient 
Fourier  transform  technique  to  obtain  the  solution  to  the  scalar  wave  equa¬ 
tion.  When  the  problem  is  discretized,  it  becomes  necessary  to  incorporate 

certain  sampling  theorems  to  ensure  an  accurate  numerical  result.  Criteria 
have  been  derived  (Ref.  2)  which  relate  the  accuracy  of  the  scalar  wave  equa¬ 
tion  solutions  for  a  desired  propagation  step  (defined  by  a  particular 
collimated  Fresnel  number)  to  the  number  of  sampling  points,  and  to  the  width 
of  the  discrete  calculation  region  relative  to  the  initial  beam  size. 

To  satisfy  these  requirements,  the  SOQ  code  user  should  describe  the 
beam  calculation  region  and  number  of  points  by  simultaneously  solving  the 
following  two  inequalities.  The  width  of  the  calculation  region  Dca^c 
required  to  model  a  propagating  collimated  beam  with  initial  diameter  2a  must 
be  chosen  to  satisfy 

1 

G  >  1  +  — r -  (3 

2  it  N  s 
c 

2 

where  G  *  DCalc^a  is  t*ie  ^arci  3  a  /LA  is  the  collimated  Fresnel 


number  for  the  propagation  step  of  length  L,  a  is  the  radius  of  the  geometric 
beam,  A  is  the  wavelength  of  the  propagating  wave,  and  £  is  the  tolerance  or 
fractional  power  spillover  outside  the  calculation  region.  For  general  SOQ 
calculations,  s  is  usually  chosen  to  be  0.02.  The  number  of  sampling  points 
N  in  each  transverse  dimension  must  be  chosen  large  enough  to  encompass  the 
spatial  frequency  bandwidth  of  the  propagating  wave.  In  Reference  2,  it  is 
shown  that 

N  >  AG  (G  ♦  1)  N 
p  —  v  c 

where  G  is  given  by  inequality  (3)  and  is  the  Fresnel  number  defined 
above.  For  GDL  SOQ  calculations,  the  value  of  *  128  is  usually  selected. 

As  an  additional  precaution,  a  raised  casine  data  window  is  employed  in  a 
user  specified  region  of  D  in  both  spatial  frequency  space  and  in  real 


2.  Sziklas,  E.A.,  and  Siegman,  A.E.,  ’’Mode  Calculations  in  Unstable  Resona¬ 
tors  With  Flowing  Saturable  Gain.  2:  Fast  Fourier  Transform  Method/’ 
Applied  Optics,  Vol  14,  pp.  1874-1899,  August  1975. 


12 


space  to  avoid  the  errors  introducted  by  aliasing  and  leakage.  Both  spatial 
frequency  and  real  space  windows  are  usually  chosen  as  ten  percent  of  the 
total  calculation  region  in  their  respective  spaces. 


Since  each  propagation  of  the  wave  nav  be  expressed  in  terms  of  the 
propagation  of  an  equivalent  collimated  propagation  step,  the  user  has  a 
means  of  determining  the  relevant  numerical  parameters  to  ensure  numerical 
accuracy  for  the  problem  being  simulated. 

4.  GENERAL  SOQ  CODE  STRUCTURE 

The  SOQ  code  can  be  viewed  as  having  several  levels  of  program  communi¬ 
cation.  Each  level  is  associated  with  executive  and/or  computational  respon¬ 
sibility.  Executive  routines  denote  those  programs  which,  based  on  user 
input,  select  and  integrate  the  appropriate  subprograms  to  effect  the  desired 
computations.  There  are  essentially  four  levels  within  this  hierarchy  which, 
macroscopically ,  describe  the  SOQ  code  structure.  The  lower  the  level  number, 
the  more  executive  is  the  routine.  The  higher  the  level  number,  the  less  exec¬ 
utive  and  more  computational  the  routine  becomes.  Care  has  been  exercised  in 
an  effort  to  provide  as  little  overlap  as  possible  between  the  levels.  The 
following  diagram  defines  the  hierarchy  which  exists  within  the  SPO  code. 


Level 

I 

SOQ 

Level 

II 

GDL 

QUAL 

LIST80 

Level 

III 
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REGAIN 

TILT 

APRTR 
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CAVITY 

OUTPUT 

MODER 

MIRROR 

STEP  SILVER 

RGRD 

TBLOOM 

THERMAL 

BLUMIT 

I  PLOT 
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Level  IV  GAINXY 

K1NET 
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SPT.AN 

DENSY 

INTERP 

POWWOW 

SIMPGG 

RAN  DU 

MIX 

POWR 

FOURT 

LINTERP 

ERF 

BLUMIT 

FUHS 

ROSN 

ERFC 

Routines  in  Level  I  are  used  primarily  to  set  ut>  the  calculations  «*t  tl 
initialization  of  a  SOQ  run  and  to  provide  flags  for  higher  level  computa¬ 
tions  and  branching.  The  program  at  Level  I  communicates  with  Levels  II 
through  Level  IV,  and  also  with  the  user  through  his  input  job  stream.  At  the 
single  Level  I  routine  "SOQ,"  the  user  decides,  after  specifying  and  con¬ 
structing  a  complex  input,  whether  far-field  (QUAL)  or  optical  train  (GDL) 
calculations  are  to  be  performed. 

Level  II  routines  are  primarily  executive  routines,  although  many  com¬ 
mutations  are  performed  in  this  level.  Subroutine  GDL  is  used  to  communicate 
with  the  user's  namelist  input,  to  direct  calls  to  Level  III  optical  train 
routines,  and  to  perform  numerical  computations.  Subroutine  QUAL,  on  the 
other  hand,  is  used  primarily  to  calculate  the  beam's  far-field  performance. 
Subroutine  LIST80  appears  at  Level  II  since  it  primarily  echoes  the  user 
input  for  a  particular  run,  and  can  only  be  accessed  through  Level  I. 

Routines  in  Level  III  are  the  principle  group  of  routines  which  a  user 
may  access  and  are  the  routines  which  make  up  the  library  of  standard  300 
optical  transmission  functions.  These  normally  are  generally  called  from 
within  GDL  and,  as  such,  receive  their  input  primarily  through  the  namelist 
input  provided  for  subroutine  GDL.  Subroutine  CAVITY  is  the  only  exception 
because  it  has  its  own  namelist  which  is  read  from  subroutine  CAVITY  the 
first  time  a  particular  cavity  element  is  called. 

Level  IV  routines  are  primarily  numerical  computation  routines  which  are 
accessed  by  Level  II  and  III  routines.  Level  IV  routines  are  used  primarily 
to  perform  numerical  calculations,  with  little  executive  responsibility. 

System  computations  based  on  namelist  input  are  primarily  effected  by 
transfer  of  data  through  named  commons,  with  subsequent  storage  of  innut  data 
in  common  to  satisfy  the  need  for  program  calculational  memory  when  iterative 
calculations  (resonator)  are  required. 
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To  meet  specific  needs,  additional  optical  element  models  may  be  con¬ 
structed.  The  user  usually  incorporates  the  new  model  in  the  Level  III 
library  of  routines.  Subroutines  called  by  the  model  are  generally  placed  in 
Level  IV,  and  modification  of  the  Level  II  program  is  generally  advised  to 
provide  input  and  complex  field  communication  to  the  new  Level  III  subpro¬ 
gram. 


3.  SYSTEM  MACROSCOPIC  FLOW  CHART 

The  system  macroscopic  flow  chart  is  the  principal  diagram  used  by  the 
SOQ  code.  The  structure  acquaints  the  user  with  the  way  information  is  passed 
from  input  to  each  subroutine  within  the  SOQ  code  along  with  a  synopsis  of 
the  method  used  to  store  information  for  subsequent  iterations.  Figure  5 
depicts  the  general  SOQ  flow  of  calculations. 


First  Level 
of 

Branching 


Figure  5.  General  SOQ  flow  of  calculations  with  memory 
arrays  associated  with  each  subroutine. 


At  the  main  program  level,  input  is  provided  which  directs  the  S00  code 
to  either  perform  beam  quality  calculations  or  resonator  and/or  optical  train 
performance  calculations.  The  parameter,  NCALL,  in  the  SOQ  main  program  input 
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argument  list  specifies  the  branch  to  be  followed.  Within  either  subroutine 
GDL  or  QUAL  the  appropriate  branches  are  taken  based  on  input  at  these  levels 
which  describe  the  input  parameters  for  the  requested  optical  transmission 
function  element. 

Since  calculations  may  be  of  an  iterative  nature,  the  SOQ  code  within 
the  GDL  subroutine  provides  an  array  of  information  which  stores  the  order 
and  name  of  each  optical  element  along  with  the  appropriate  input  that  char¬ 
acterizes  that  element.  Correspondence  between  the  equivalent  optical  train 
which  the  optical  analyst  has  described  is  retained  within  the  memory  system 
of  the  SOQ  code.  Figure  5  depicts  the  structure  of  the  SOQ  code  with  the 
memory  arrays  associated  with  each  subroutine. 

For  each  GDL  subroutine  initially  called,  the  model  input  is  retained 
through  the  IGDL  array  shown  in  Figure  6.  For  each  model,  the  ABC  (P^,  Q^, 

R^)  array  retains  the  input  parameter  values.  Iterative  solution  of  the  reso¬ 
nator  problems  uses  the  values  stored  in  the  IGDL,  ABC,  or  similar  arrays  to 
direct  the  optical  train  through  the  same  path.  Noniterative  calculations  are 
also  stored  in  the  IGDL  and  ABC  arrays,  although  the  arrays  are  not  used  to 
redirect  calculations. 


Figure  6.  Subroutine  GDL. 


The  notable  exception  to  this  philosophy  is  subroutine  CAVITY  with  its 
own  input  and  memory  system  (Fig.  7).  The  subroutine  CAVITY  arrays  are 
structured  according  to  the  cavity  input,  and  are  a  function  of  the  number  of 
cavities  used. 


SOQ 


Figure  7.  Subroutine  CAVITY  with  its  own  input  and  memory 
system. 

The  overall  SOQ  macroscopic  flow  chart  can  now  be  constructed.  The  over¬ 
all  flow  chart,  shown  in  Figure  3,  delineates  the  main  program  and  the  sub¬ 
routine  called  from  each  level.  Basically,  two  branches  ("QUAL”  and  "GDI/') 
can  be  followed.  Within  each  branch  the  precise  calculations  and  input  are 
read  to  define  the  physical  problem  being  simulated. 

6.  OVERVIEW  OF  SOQ  CODE  OPERATION 

A  previously  defined  optical  train  with  known  relevant  physical  para¬ 
meters,  the  analyst  may  identify  the  SOQ  routines  necessary  to  simulate  the 
optical  train  or  configuration  of  interest.  With  the  relevant  SOQ  model (s), 
the  analyst  may  use  the  detailed  description  of  the  model  contained  herein  to 
correlate  his  physical  input  parameters  with  those  required  to  run  the  SOQ 
code . 
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The  input  is  then  coded  with  SOQ  namelist  input  descriptions  to  direct 
the  calculations  to  the  appropriate  subprograms.  The  following  demonstrates 
the  procedure  that  may  be  used  to  execute  the  SOQ  code.  Actual  test  cases 
specifically  illustrating  SOQ  code  setup  are  included  in  Section  V  and  are 
maintained  on  the  AFWL-CDC  Cyber  176  computer. 
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I 


Initialize  SOQ  parameters 
and  load  input 
distribution 

II  Transfer  to  "GDL"  for 
optical  train/ resonator 
calculations  or  calculate 
beam  quality  through  "QUAL" 

III  In  GDL  read  namelist  param¬ 
eters  to  pick  model 

IV  When  model  has  been  de¬ 
fined,  GDL  branches  to 
the  appropriate  namelist 
for  the  chosen  model 

V  Repeat  input  calls  to 
initialize  additional 
models 

VI  Branch  out  of  GDL  and  stop, 
or  start  a  new  set  of 
optical  model  calculations. 


Namelist  START  activates 
code  "SOQ" 

Parameter  NCALL  =  2  in 
START  namelist  activates 
"GDL"  NCALL  =  3  activates 
"QUAL" 

Namelist  CONTROL,  IFLOW  * 
IVWX 

IVWX  =  Chosen  model  number 

Namelist  XYIPQ,  where 
XYZPQ  is  associated  with  an 
appropriate  chosen  model , 
such  as  an  aperture  or 
mirror 

Repeat  steps  II  through  IV 

Namelist  CONTROL,  IFLOW  = 

9,  returns  to  program 
SOQ  from  GDL 


The  actual  job  setup  necessary  to  achieve  this  is  described  schematically 
in  Figure  9. 

A  more  detailed  analysis  with  examples  of  the  input  and  permanent  (perm) 
file  stream  necessary  to  execute  the  "SOQM  code  follows  this  illustration. 

7.  USER  ACCESS  AND  PERMANENT  FILES  REQUIRED 


Input : 

The  SOQ  code  is  configured  so  that  all  input  parameters  are  input 
through  standard  namelists.  To  find  the  specific  namelist  required,  go  to 
Section  III  where  inputs  for  the  chosen  subroutines  are  identified.  Section 
III  describes  the  basic  procedure  used  and  namelist  called  to  initiate  the 
optical  calculations.  This  is  best  illustrated  in  Figure  10. 


Figure  9 


Job  setup  necessary  for  providing  subprograms 
of  interest. 
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Figure  10.  Section  III  describing  the  basic  procedure  used 
and  the  requisite  namelist  called  to  initiate 
optical  calculations. 

If  the  user  desires  quality  calculations  in  the  branch  to  QUAL.  the  code 
looks  for  input  provided  by  namelist  QLOT.  If  the  user  chooses  the  GDL  cal¬ 
culation  flowpath,  then  the  code  expects  its  next  level  of  input  from  name- 
list  CONTRL  associated  with  GDL.  After  close  inspection  of  subroutine  GDL, 
the  user  will  determine  the  appropriate  namelist  to  specify  his  optical 
train.  As  each  new  model  is  applied  in  the  optical  train,  the  user  executes  a 
new  }  CONTRL  input  followed  by  the  appropriate  model/transmission  function 
input.  The  following  computer  printout  is  included  to  demonstrate  the  calling 
sequence.  Without  particular  regard  for  the  parameters  used,  the  user  is 
directed  to  the  sequence  of  inout  and  the  comments  appended.  In  this  set  of 
input,  Item  I  input  is  used  to  initialize  the  SOQ  main  program  and  branch  to 
subprogram  GDL.  Items  II  through  VII  control  branches  within  GDL.  Item  VIII 
branches  out  of  GDL  back  to  the  main  program  and  Item  IX  branches  in  SOQ  to 
terminate  the  SOQ  run.  All  SOQ  input  is  handled  in  the  same  manner,  whether 
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8.  PERMANENT  FILE  REQUIREMENTS  AND  JOB  EXECUTION 

During  execution  the  SOQ  computer  program  may  use  several  random  access, 
on-line,  permanent  files.  These  files,  called  tapes  but  actually  disk  files, 
are  used  to  access  the  data  base  or  write  to  the  data  base.  Permanent  files 
or  scratch  files  may  be  created  by  storing  complex  field  information  repre¬ 
senting  the  propagating  wave,  the  complex  gain  distribution  in  cavities,  or 
far  fields  distribution  information.  Internally,  the  code  also  stores  complex 
field  information  on  files  to  protect  it  during  some  operations.  This  discus¬ 
sion  is  included  to  acquaint  the  user  with  the  hierarchy  of  files  which  are 
reserved  and  freed  during  the  SOQ  code  execution. 

Reference  is  made  to  the  program  definition  in  Program  SOQ  which  defines 
the  tapes  used  by  the  program. 

Additionally,  the  user  must  attach,  as  a  minimum,  the  SOQ  code  either  as 
a  source  file  or  as  a  binary  file  to  provide  a  compilable  and/or  exectuable 
module  for  the  computer. 

The  code  currently  defines  tapes  1  through  31  in  the  program  card.  The 
following  delineates  the  file  correspondence  for  normal  SOQ  calculations. 

TAPE _ 

1  Used  automatically  by  QUA1 

2  Open 

3  Open 

4  Open 

5  Reserved  for  input  file 

6  Reserved  for  output  file 

7  Reserved  for  scratch  storage 

8  Reserved  for  cutout  of  feedback  resonator 

field 

9  Reserved  for  complex  field  and  associated 
coordinates 

10  Used  automatically  by  CAVITY 

11  Reserved  for  cavity  calculations 

12  Reserved  for  cavity  calculations 

13  Reserved  for  cavity  calculations 

15  Reserved  for  cavity  calculations 
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TAPE 


16  Reserved  for  cavity  calculations 

IT  Reserved  for  cavity  calculations 

13  Reserved  for  cavity  calculations 

19  Reserved  for  cavity  calculations 

20  Reserved  for  cavity  calculations 

21  Used  automatically  by  CAVITY 

22  Reserved  for  cavity  calculations 

23  Reserved  for  cavity  calculations 

24  Reserved  for  cavity  calculations 

25  Reserved  for  cavity  calculations 

26  Reserved  for  cavity  calculations 

27  Reserved  for  cavity  calculations 

28  Reserved  for  cavity  calculations 

29  Reserved  for  cavity  calculations 

30  Reserved  for  digital  density  field 

31  Reserved  for  digital  density  field 

The  user  may  increase  the  allocation  of  logical  files,  but  must  be  aware 
of  the  usual  file  correspondence  that  is  currently  automatically  allocated  in 
the  program. 

The  attached  day  file  for  a  normal  SOQ  resonator  run  establishes  the 
usual  file  correspondence  and  delineates  the  way  the  user  obtains  an  execut¬ 
able  module  under  the  NOS/BE  Cyber  176  system.  The  following  computer  print¬ 
out  is  an  example. 

9.  SOQ  OUTPUT  FORMAT 

When  properly  executed,  the  SOQ  code  will  calculate  relevant  system  per¬ 
formance  either  in  the  near  field  (Fresnel)  or  far  field  (Fraunhofer)  sense. 
Each  of  the  SOQ  transmission  functions  which  modifies  the  complex  field  is 
programed  to  display  or  echo  the  relevant  input  parameters,  the  calculated 
relevant  system  parameters,  and  the  amount  of  energy  in  the  beam.  The  rou¬ 
tines  generally  calculate  integrated  energy  loss,  if  applicable,  and  the 
location  of  maximum  or  peak  flux  loads.  The  description  of  each  subroutine, 
along  with  its  listing  orovides  insight  into  the  actual  output  values.  The 
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system  of  test  cases  provided  on  the  AFWL/CDC  computer  allows  the  user  to 
exercise  the  SOQ  routines  and  to  evaluate  the  output  parameters  with  the 
individual  routines. 
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By  providing  performance  information  at  each  optical  element,  the  user 
may  determine  those  specific  areas  at  each  element  which  violate  his  system 
constraints.  The  SOQ  code  also  provides  far-field  performance  calculations, 
through  subprogram  QUAL.  The  output  of  these  calculations  allows  the  user  to 
assess  the  far-field  performance  of  the  complex  field  distribution.  The  test 
cases  provided  in  Section  V  demonstrate  the  far-field  calculation  output. 

10.  SOQ  ROUTINES 

SOQ  routines  are  presented  below  to  establish  correlations  between  the 
high  energy  laser  modeling  task  and  the  names  of  the  specific  SOQ  routines 
within  specifically  identifiable  areas.  The  routines  may  be  called  up  to  per¬ 
form,  either  in  whole  or  in  part,  the  desired  modeling  function.  A  more 
detailed  description  of  each  model  or  element  is  presented  in  subsection  K, 
the  Director  of  SOQ  Models. 


GENERAL  AREAS: 

SUBROUTINES 

Laser  Cavity  Calculations: 

Kinetics  of  lasing  media 

CAVITY,  KINET,  GAINXY, 

SIMPGG,  MIX 

REGAIN 

Density  I homogeneities 

DENSY,  AEROW 

Beam  Propagation 

STEP,  GDL,  CAVITY 

FOURT 

Far  Field  Calculation 

QUAL,  PLTOT,  POWWOW 

Tilt,  CENBAR 

Numerical  Interpolation 

INTERP,  ROSN ,  R0SN6 

LINTERP,  RGRD 

Thermally  Induced  Distortions 

BLUMIT,  TBLOOM 

THERML,  FUHS 

Optical  Elements 

Resonator  Iteration  S  Mode  Proper¬ 

MIRROR,  APRTR,  SLIVER 

SPIDER 

ties  within  GDL 

N10DER 

General:  Write  or  echo  input  to 

obtain  printer  plots  of  complex 

distribution 

LIST80 ,  OUTPUT,  IPLOT 
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11.  DIRECTORY  OF  SOQ  MODELS 


SOQ 


GDL 


QUAL 


AEROW 


Executive  program  used  to  initialize  field  distribution,  call 
optical  train  calculations  through  subroutine  GDL,  or  compute  far 
field  performance  through  calls  to  subroutine  QUAL. 

Executive  and  computational  program  used  to  direct  optical  train 
and  resonator  calculations.  Subroutine  GDL  may  perform  any  of  the 
following  tasks: 

1 .  Read  input 

2.  Apply  loaded  cavity  calculations 

3.  Apply  mirror  model 

4.  Apply  propagation  algorithm 

5.  Apply  aperture 

6.  Apply  thermal  blooming 

7.  Determine  convergence 

S.  Plot  field  distributions 

9.  Branch  to  program  SOQ 

10.  Load  or  read  field  from  disk 

11.  Apply  aerodynamic  window  model 

12.  Orthogonal ize  resonator  modes 

13.  Scale  entire  complex  field  by  arbitrary  magnifica¬ 
tion 

14.  Flip  field  about  chosen  axis 

15.  Apply  sinusoidal  density  variations 

16.  Regrid  the  field  to  different  sized  DCALC 

17.  Write  field  to  punched  cards 

18.  Apply  mirror  thermal  boundary  layer  model 

19.  Apply  telescope  spider 

20.  Perform  Fox-Li  iterative  resonator  calculations. 

Executive  and  computational  routine  for  the  calculation  of  far 
field  complex  distribution.  Subroutine  QUAL  may  be  used  to: 

1.  Identify  the  tilt  and  defocus  components  in  optical  field 

2.  Computer  the  field  centroid 

3.  Determine  the  far  field  irradiance  distribution 
Apply  the  aerodynamic  window  root-mean-square  (RMS)  phase 


RANDU 

APRTR 

BLUMIT 

CAVITY 

CENBAR 

DENSY 

FOURT 

FUHS 

GAINXY 

INTERP 

IPLOT 

KINET 

MIRROR 

MIX 


W 


Calculate  uniformly  distributed  random  numbers  on  the 
interval  (0,1). 

Apply  various  types  of  apertures 

1.  Circular 

2.  Rectangular 
5.  Annular 

Apply  and  calculate  intra-  or  intercavity  thermally-induced 
phase  distortion  due  to  interaction  of  the  complex  field 
with  the  cavity  gases. 

Apply  the  effect  of  nonlinear  saturable  GDL  gain  and  phase 
disturbance  to  a  complex  field  and  propagate 

Calculate  the  beam  centroid. 

Calculate  and  load  the  density-induced  phase  distortion 
for  cavity  calculations. 

Calculate  forward  or  backward  Fourier  Transform  of  a 
specified  two-dimensional  input  distribution  via  Fast 
Fourier  Transform  techniques. 

Calculate  and  apply  the  effects  of  lower  laser  level-to- 
ground-state  heat  addition  perturbations  to  cavity  den¬ 
sity  in  the  supersonic  flow  field. 

Calculate  small  signal  or  loaded  gain  for  cavity. 

Perform  linear  interpolation  of  an  arbitrary  complex 
field  distribution. 

Create  line  printer  plots  of  the  complex  field  including 
plots  of  the  intensity  and  phase. 

Calculate  kinetics  and  small  signal  gain  for  use  in  the 
routine  GAINXY. 

Apply  the  effects  of  mirrors  to  a  propagating  field,  in¬ 
cluding  nonlinear  interaction  of  the  complex  field  and 
the  induced  mirror  surface  response. 

Calculate  relaxation  and  pumping  rrtes  for  KINET, 
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MODER 


Apply  mode  orthogonalization  capability  via  Gramm-Schmidt 
procedure  to  calculate  higher  order,  bare,  or  loaded 
cavity  modes. 


OUTPUT 

PLTOT 

POWWOW 

REGAIN 

RGRD 

ROSN 

LINTERP 

ROSN  6 

SIMPGG 

SILVER 

SPIDER 

SPTAN 

STEP 


Generate  amplitude  C or  intensity)  and  phase  printer  plots  at  three 
slices  through  the  field. 

Calculate  and  produce  printer  plot  of  far  field  distribu¬ 
tion  and  calculate,  via  POWWOW,  the  far  field  integrated 
irradiance  within  a  specified  far  field  radius. 

Calculate  intercepted  power  (integrated  irradiance)  in 
far  field  coordinates. 

Calculate  and  apply  the  effect  of  the  complex  field  on 
the  nonlinear  gain  medium  during  resonator  iterative 
calculations. 

Regrid  the  current  SQQ  complex  field  to  a  user-defined 
size. 

Perform  cubic  spline  interpolation  of  interferometric 
density  data  stored  on  peripheral  device. 

Perform  linear  interpolation  of  density  field  dati 
parameterized  by  aerodynamic  sidewall  projection. 

Bivariate  interpolation  of  the  spline  fit  data  using 
cubic  splines.  Allow  inclusion  of  the  cavity  density 
field  from  ‘direct  interferogram  data  reduction. 

Calculate  loaded  gain  by  3-level  closed  form  solution  for 
a  N2-C0?  GDL  system. 

Apply  an  annular  aperture  to  the  complex  field. 

Apply  a  multi-element  beam  obscuration  as  found  in  a 
telescope  spider. 

Obtain  a  tangent  function  on  the  interval  (0,  2r) . 

Propagate  the  complex  field,  via  scalar  wave  theory, 
using  a  variable  or  constant  mesh.  Fourier  Transform 
is  used  to  solve  the  scalar  wave  equation. 
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TBLOOM 


Apply  and  calculate  nonlinear  thermally-induced  phase 
distortion  on  the  complex  field  due  to  interaction  of  the 
complex  field  with  the  medium  through  which  it  is 
propagating. 

THERMt  Apply  and  calculate  effects  of  thermal ly- induced  boundary 
layers  on  the  complex  field. 

TILT  Calculate  the  complex  field  tilt  and  sphere. 

ERF,  ERFC  Evaluate  the  error  function  and  complementary  error 
function. 

12.  FUNDAMENTAL  COMPUTATIONAL  STRUCTURE  OF  THE  SOQ  CODE 

The  purpose  of  this  section  is  to  define  the  correspondence  between  the 
complex  numerical  field,  CU,  used  in  the  SOQ  code  and  the  real  Cartesian 
coordinate  system.  Throughout  the  SOQ  code,  the  fundamental  approach  to  cal¬ 
culating  the  complex  field  has  been  to  represent  the  optical  field  interac¬ 
tion  as  a  complex  transmission  function.  The  complex  field,  then,  is  gener¬ 
ated  and  modified  at  each  optical  element/subelement  in  the  following  format’ 

Or  (x,y)  =  CU(x,y)  *  T  (x,y,c,  CU(x,y),  P)  (5 


where  x,v  represent  transverse  spatial  coordinates  (cm),  z  represents  coor¬ 
dinates  along  the  optical  axis  (cm),  CU(x,y)  represents  the  complex  field  in 
Cartesian  (a(x,y),  b(x,y))  form,  and  represents  a  vector  of  scalar  param¬ 
eters  specific  to  the  transmission  function  model,  T,  being  applied.  T 
represents  the  complex  transmission  function. 

It  is  assumed  in  the  SOO  computer  code  that  each  operation  on  the  field, 
CU,  can  be  represented  as  a  complex  transmission  function  or  as  a  product  of 
transmission  functions.  Each  optical  operation  or  "interaction"  of  the  prop¬ 
agating  field  with  the  surroundings  modifies  the  transmission  function  frame¬ 
work. 

Within  the  SOQ  code,  the  primary  carrier  of  this  field  information  is 
the  complex  array  CU,  stored  in  common  /MELT/.  The  correspondence  between 
ohysical  coordinates  and  the  CU  array  indices  is  illustrated  in  Figure  11. 
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+DCALC.  DCALC 


Correspondence  between  physical  coordinates  and  Cll  array  iiulices. 


The  calculation  region  is  defined  to  be  a  DCALC  X  DCALC  square  sub¬ 
divided  by  NPTS  number  of  points  to  define  the  soatial  coordinates,  where 
MPTS  *  2K  i.e. . 

AX  =  DCALC/NPTS  (6) 

DCALC  AX 
1  2  2 


X2  =  *1  +  AX 

Xn  *  Xi  «■  (n-1)  AX 


(71 


Since  the  mesh  is  square,  no  additional  storage  for  y  is  defined  because  the 

x  and  y  coordinates  may  be  represented  by  the  same  array.  Additionally,  the 

2 

CU  array  is  defined  as  a  singly  dimensioned  string  whose  length  is  (NPTS) ”  = 
MOB.  The  correspondences  of  CU  with  the  real  spatial  coordinates  are: 


01(11 

CU(2) 

CU(3) 

CU  (MPTS1 
CU  (NPTS+1) 
CU  fNPTS+2) 
etc. 


X(l),  Y(l) 
X(2) ,  Y  ( 1) 
X(3),  Y(l) 

X(NPTS),  Y(l) 
X(l),  Y(2) 
X(2) ,  Y(2) 


(X(I)  =  Y(D) 


The  CU  array  and  the  X  array  are  stored  in  common  /MELT/.  Other  arrays 
stored  in  /MELT/  are: 


CFIL 

Dummy  work  array 

WL 

Radiation  wavelength 

NPTS 

Number  of  points  in  x  direction 

DRX 

Integrated  displacement  of  optical 

line  in  x  direction. 

axis 

from  beam  center- 

DRY 

Integrated  displacement  of  optical 
line  in  y  direction. 

axis 

from  beam  center- 
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Principally,  common  /!!ELT/  is  the  carrier  of  the  comnlex  field  informa¬ 
tion  at  any  location  of  station  in  the  optical  train.  The  correspondences 
between  the  Cartesian  (x,y,3)  coordinate  system  and  the  indices  of  the  CU  and 
X  arrays  are  maintained  throughout  the  SOQ  code  as  a  standard  field  orierta- 
tion  scheme,  within  which  each  optical  transmission  function  is  aoolied. 


